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Abstract. We study systems of Brownian particles on the real line, which interact 
by splitting the local times of collisions among themselves in an asymmetric manner. 
We prove the strong existence and uniqueness of such processes and identify them 
with the collections of ordered processes in a Brownian particle system, in which 
the drift coefficients, the diffusion coefficients, and the collision local times for the 
individual particles are assigned according to their ranks. These Brownian systems 
can be viewed as generalizations of those arising in first-order models for equity 
markets in the context of stochastic portfolio theory, and are able to correct for 
several shortcomings of such models while being equally amenable to computations. 
We also show that, in addition to being of interest in their own right, such systems 
of Brownian particles arise as universal scaling limits of systems of jump processes 
on the integer lattice with local interactions. A key step in the proof is the analysis 
of a generalization of Skorokhod maps which include 'local times' at the intersection 
of faces of the nonnegative orthant. The result extends the convergence of TASEP 
to its continuous analogue. Finally, wc identify those among the Brownian particle 
systems which have a probabilistic structure of determinantal type. 



I. Introduction 

Systems of Brownian particles with various types of interactions have been widely 
studied. More recently, Brownian particles with electrostatic repulsion (Dyson's 
Brownian motion, see for instance section 4.3 in [1]) play a central role in the un- 
derstanding of the universality properties of large Hermitian random matrices with 
independent entries (Wigner matrices, see the survey [7j and the references there). 
In addition, systems of Brownian particles interacting through their ranks, which 
were originally introduced in the context of the piecewise-linear filtering problem and 
in the resulting study of diffusions with piecewise-constant characteristics (see [1]), 
have been of great importance in the study of large equity markets within stochastic 
portfolio theory (see [8], [ID]). Finally, systems of Brownian particles interacting by 
stickiness have been recently introduced as continuous analogues of a certain ran- 
dom evolution for the distribution of mass on the integer lattice (see [16] and the 
references there). 



Date: October 2, 2012. 

2000 Mathematics Subject Classification. 60K35, 60H10, 91B26. 

Key words and phrases. Determinantal processes, Interacting particle systems, Invariance prin- 
ciples, Reflected Brownian motions, Skorokhod maps, Stochastic Portfolio Theory, Strong solutions 
of stochastic differential equations, Triple collisions. 

Research partially supported by NSF grants DMS-09-05754 (I. Karatzas) and DMS-10-07563 (S. 
Pal). 

1 



2 



IOANNIS KARATZAS, SOUMIK PAL, AND MYKHAYLO SHKOLNIKOV 



As a starting point towards the formulation of our setup, we recall from section 
13.1 in [10] that the ordered processes in a system of Brownian particles interacting 
through their ranks are given by independent Brownian motions, which have con- 
stant drift and diffusion coefficients and collide in a symmetric fashion; that is, each 
collision local time is split equally between the two colliding particles. In contrast, 
the particles in the process introduced by Warren in [25J evolve as independent 
Brownian motions which have constant drift and diffusion coefficients and collide in 
a totally asymmetric manner; that is, the collision local time is assigned entirely to 
one of the two colliding particles. 



1.1. The Setup. The dynamics of ( 1.1 ) below give the formal description of a system 
of ordered Brownian particles on the line, which move as independent Brownian 
motions with constant drift and diffusion coefficients and collide asymmetrically; 
that is, the collision local times are apportioned unequally, in a manner that depends 
on the ranks of the particles involved in the collisions. 

Consider a continuous, n— dimensional semimartingale (-Ri(-), R2 (•))■■•> Rn(')) 
with values in the Weyl chamber W n = {(r±, r2, . . . , r n ) : 00 > 7*1 > r-i > • • • > r n > 
—00} , and with dynamics of the form 

(1.1) R k (t) = R k ^) + b k t + a k (3 k {t) + q- ^ k ' k+l \t)-qt ^ k - 1 ' k \t), < t < 00 

for k = 1, 2, . . . , n . Here the drifts b\, 62, • • • , b n are given real numbers; the dis- 
persions o"i, o"2, . . . , o~ n are given positive real numbers; the collision parameters 
qf , qf,...,q^ are given positive real numbers satisfying 

(1.2) q~ +q+ +1 = 1, fc = 1,2,..., Ti-1; 

and the processes /3 1 (-), f3 2 (•)>••• , /3 n (-) are independent standard Brownian mo- 
tions. On the other hand, for each k = 1,2, . . . ,n — 1 the process 

(1.3) A (fc ' fc+1) (-) = L Rk - Rk +^-;0) 

is the right-sided local time accumulated at the origin by the nonnegative semimartin- 
gale R k (-)-R k+1 (-) ; we set A^-) = A(™."+ 1 )(.) =0. The "regulating" role of these 



local times in (1.1) is to make sure the resulting process (-Ri(-), ^(O; ■ ■ > Rn(')) =: 
R(-) takes values in the wedge W™ at all times. This process R(-) can thus be 
regarded as Brownian motion with reflection on the faces of the polyhedral domain 
W n , in the sense of Harrison & Williams (1987) (cf. [T3]. [2B]k 



As we show in section 3, Brownian particle systems of the type (1.1) are not only 
of interest in their own right, but also arise as universal scaling limits for systems of 
jump processes on the integer lattice with local interactions. Consider a system of n 
particles on the integer lattice, moving according to (possibly asymmetric) continuous 
time simple jump processes, which are independent as long as the particles are located 
at n different sites. When two or more particles land at the same site (we will refer to 
such events as "collisions"), the jump rates of the particles change, and in a manner 
that preserves the order of the particles. As we explain in section 3, such particle 
systems converge under a diffusive rescaling of time and space to the solution of a 



stochastic equation of the type (1.1). The proof of this result necessitates a detailed 
study of Skorokhod maps that transform noise to processes constrained to stay in 
the nonnegative orthant, but might involve 'local time' push from the intersection of 
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multiple faces of the orthant. Along the way, we generalize the boundary property of 
reflected Brownian motion established in [22] (see Lemma[l]below) and the invariance 
principle for reflected Brownian motion of (27] (see Proposition M below). 



We now discuss how the processes in (TTTj) can be seen as describing the order 
statistics in Brownian particle systems, in which the particles are allowed to ex- 
change their ranks. As we explain below, the latter can be used as models for the 
logarithmic capitalizations in large equity markets, and generalize the so-called "first 
order models" of stochastic portfolio theory. Consider an n— dimensional process 
(Xx(-), X 2 (-), . . . , A„(-)) that satisfies the system of stochastic differential equations 



dXi(t) = J2 1 {x t (t)=R^(t)} b kdt + J2 1 {xdt)=R^(t)} cr kdW i (t) 

k=l k=l 

n 

(1-4) + £ WHHfW) " (V2)) dA^ +1 >(t) 

k=l 

n 

-EW (g fc + -(l/2))dA^^)(t). 

k=l 

Here Wi('), W 2 (-), ■ ■ ■ , W n (-) are independent standard Brownian motions; the "de- 
scending order statistics" 

(1.5) max X 3 {.) =: R* (■) > R* (•) > > := min Xj(-) 

are the ranked processes corresponding to the particle system X 1 (-),X 2 (-), . . . , X n (-) 
with lexicographic resolution of ties; whereas 

(1.6) A^(-) = L R k~ R f(- ; 0) , t > k + 1 

denotes the local time accumulated at the origin by the nonnegative semimartingale 

R^-)-Rf(-). 



Assume now that a weak solution to the system of (1.4) exists and satisfies 

(1.7) Ceb ({t > | 3 1 < % < j < n : X^t) = X^t)}) = 

almost surely, where Ceb denotes the Lebesgue measure on [0, oo). Then with the 
notation J\ff.(t) = \{i : Xi(t) = R*(t)}\ for the number of particles occupying the 
k-th. rank at time t , the Banner & Ghomrasni (2008) formula 
(1.8) 

n / n k—l N 

ditf(*) = Eh^t ) =x m dX i (t) + - m £ dA^(t)-]TdA^)(t) 
i=i ' \e=k+i t=i , 

(cf. Theorem 2.3 in [3]) shows that the process of spacings 

(1-9) Z x (.) := {R?(.)-RX(.),R*(.)-RX(.) r .. , R^(-) - R* (■)) , 

when away from the boundary of the nonegative orthant (R + ) n_1 , moves according 
to the multidimensional process ((&i — b 2 )t + (j\ (3f(t) — a 2 P 2 (t), . . . , (b n ^\ — b n )t + 
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°n-i Pn-i(t) ~ a nf3n(t)) , t > 0. Here the processes 

n „ . 

(1.10) 0* (•) := / W)=Jff(t)} dWi(t) , fc = 1, 2, . . . 



i=l 



are independent standard Brownian motions, by virtue of the P. Levy theorem (see 
section 3 in [2 J for a very similar derivation). But on the strength of Lemma [I] below, 
which generalizes the boundary property of reflected Brownian motion established by 
Reiman & Williams (1988) in [22], the triple- or higher-order collision local times 
A( fc '^)(.) vanish for all £ > k + 2 , so the expression in the Banner-Ghomrasni 
(2008) formula (|T§ simplifies to 



E 1/cv 



m dX.it) + i dA^ +1 >(t) - I dA^^t) 



i=i 

n 



{Hf(t)=jfiW} ^ A {Jfi(*)=«f (*)} 



1=1 



+ E 



=1 

n 



n 

b t dt + 1 {x l (t)=R?{t)} °l dWi(t)j 
i=i 



(t)=Xi(t)} 



EV^fC*)} fe--(V2))dA^ +1 )(0 



£=1 



E ^ w=^w} E W)=Hfw} fe + - (V2)) dA^(o 

i=l £=1 



+ I d A(M+l)( t )_^ dA (fc-U0( t ) 



= 1, . . . , n . 

Evaluating the sums of this expression over % and i and recalling the notation of 



(1.10), we get the dynamics (1.1) for the descending order statistics of (1.5). 



1.2. Interpretation. We shall think of the processes Xi(-), X 2 (-), . . . , X n {-) as rep- 
resenting the positions of a finite collection of Brownian particles, whose drift and 
dispersion coefficients are assigned according to the ranks occupied by the parti- 



cles when ordered as in (1.5). When the particles collide, they interact asymmetri- 



cally with their nearest neighbors through the collision local times at the origin, in 



the specific manner of (1.4) and with the notation of (1.6). We shall speak of the 



W n — valued semimartingale (Xi(-), . . . ,X n (-)) as the process of "names" (positions 
of individual particles), and of the components of the W- valued semimartingale 
(Ri(-), ■ ■ ■ ,Rn(')) as the associated "ranked processes" (descending order statistics). 

If we denote by t t (i) the rank occupied by particle i at time t, we can write 



(1.11) 



the system of equations (1.4) in the very informal, yet suggestive and slightly more 
compact form 

dX t (t) = b xt{{) dt + a rtW dW t (t) + (q-y - (1/2)) d A^^ +1 \t) 

- (g+ - (1/2)) dA<*»-i,*<0)(t) . 

To wit: at any given time t , every particle i gets assigned drift and dispersion 
parameters according to its current rank x t {i) , and feels an upward (respectively, 
downward) local-time-like pressure, or "drag", when colliding with the particle right 
below it (respectively, right above it) in proportion to Q^rq — (1/2) • 
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1.3. Some Special Cases. It is instructive to compare the system (1.4) in the 
two-dimensional case n = 2 with the systems of equations studied by Fernholz, 
Ichiba & Karatzas (2012) in [9] (see the systems of equations (1.2), (1.3) and 
(4.13), (4.14) of that paper). As one can see by comparing the coefficients of the 



local time terms, in the case n = 2 the system (1.4) is a special case of the system 
of equations (4.13)-(4.14) in [9], namely 

dXi(t) = l {XlW > X2W} (6 1 dt+a 1 d^ 1 (t))+l {XlW<X2W} (6 2 dt+a 2 d^ 2 (t))+/tdLl Xl - X2 l(t) 
dX 2 (t) = l {Xl(i)<X2(i)} (6 1 dt+a 1 d^ 1 (t))+l {XlW > X2W} (6 2 dt+(T 2 diy 2 (t))+fi:dLl Xl ^ 2 l(t) 

with k = gf — (1/2) = (1/2) — q 2 . The inequalities in the indicators reflect the 
convention on lexicographic resolution of ties we referred to earlier. 

In particular, with T(-) = Xt(-) - X 2 (-) , S(-) = Xt(-) + X 2 {-) , A x = h - b 2 and 
A 2 = b\ + b 2 , we have -Ri(-) — -R2(') = |T(-)| > an d the processes 

W(t) = T(t) + Ai / sgn(T(s)) ds , V{t) = E(t) - A 2 t - 2 « L |T| (£) , < t < oo 
Jo 

are now Brownian motions with diffusion coefficients \/ o\ + o\ and with covariation 
(W,V}(-) = (aj - al) / 8gn(T(t))dt. 

Let us also remark that the solution to the above two-dimensional system can be 
realized as the solution of the system of equations (1.2)-(1.3) in [9J with 

^ ' ' K Ql 2 2 ?2 2 2 2 2' 

and with rji , r] 2 , Ci, C2 the parameters introduced in [5]. In particular, the condition 
(1.5) in [9], which is necessary and sufficient for the well-posedness of the system of 
equations (1.2)-(1.3) in [9], is here fulfilled. 

• Now, suppose that we have 

b 1 = ... = b n = 0, ax = . . . = a n = 1 and g]~ = . . . = q~ = , qf = . . . = g+ = 1 . 

In this case one can think of the particles as having masses that decrease from right 
to left, so that the mass of each particle is negligible compared to the mass of its right 
neighbor. To wit, whenever a particle collides with its right neighbor, it is reflected off 
this considerably "heavier" particle. In this situation, the vector of ranked processes 
R(-) is given by the continuous version of the Totally Asymmetric Simple Exclusion 
Process (TASEP). We refer the reader to section 4 in [22] for some of the properties 
of this process, and to [UJ for its appearance as the scaling limit of TASEP. 

• Let us also note that the special case 



1 



(1-13) =2"' fc = 1 ' 2 



, n . 



in which all local times disappear from (1.4), (1.11) and get equal weights in (1.1), 



and all collisions of ranked particles are symmetric, was studied in detail by Ichiba, 
Karatzas & Shkolnikov in in this case, individual particles collide with each 
other without feeling any local time drag from their nearest neighbors. 
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1.4. Outline. The rest of the paper is organized as follows. In section 2.1 we show 



the strong existence and uniqueness of the solution to the system of equations (1.1). 



Subsequently, in section 2.2 we prove the appropriate generalization of the bound- 
ary property of reflected Brownian motion established in [22] to processes involving 
local times accumulated on lower-dimensional faces of the boundary; then, we study 
the attainability of lower-dimensional parts of the boundary by such processes in 
section 2.3. These results are combined in section 2.4 to prove the strong existence 



and uniqueness of the solution to the system of equations (1.4), under the explicit 



conditions on non-attainability of lower-dimensional boundaries derived in section 
2.3. The latter results extend those of [I], [19] and [23] (within our setup, the first 
two only address the case (1.13), whereas the last one only addresses the case of 



rank-independent parameters). Then, in section 2.5 we view the solution to (1.4) as 



a model for the logarithmic capitalizations in large equity markets and analyze the 
resulting capital distributions. 



Section 3 is devoted to the identification of the processes in (1.1) as universal 



scaling limits of systems of one-dimensional continuous time jump processes with 
local interactions. It includes an extension of the invariance principle of [27] in 
section 3.3 (see Proposition [9]). Finally, in section 4 we characterize the sets of 
collision parameters, for which the solution to (1.1) has a probabilistic structure 
of determinantal type. These results generalize those in [23], on Brownian particle 
systems with totally asymmetric collisions. 



2. Analysis of the continuous process 

2.1. Ranks. We start with the construction of the vector R(-) = (i?i(-), . . . , R n (-)) 
of ranked processes. We note first that, due to the positivity of the coefficients 
qf, qf, . . . , q^, there exist positive constants ci, C2, . . . , c„ such that Y^i=x c iRi{t)) t > 
is a Brownian motion with drift, that is, the contribution of the local times to 
its dynamics vanishes. This observation allows us to construct the ranked processes 
Ri(-) , R 2 (-) , ■ ■ ■ ,R n ( 4 ) using the following procedure: first, we define the auxiliary 
Brownian motion R(-) = i^G), • • • > -^n(')) 5 f° r which 

(n \ / n \ n 

k=l J \k=l J k=l 

and 

d(R k (t) - R k+1 (t)) = (b k - b k+1 ) dt + a k d/3 k (t) - a k+1 d/3 k+1 (t) 
for k = 1, 2, . . . , n — 1 . Next, we introduce the process 

n 

Y(-) := (Y,CkRk(-),Rl(-)-R2(-),...,Rn-l(-)-Rn(-)) 

fe=l 

and apply the Harrison- Reiman (1981) version of the Skorokhod reflection map 
ijr-ff.R f or orthant (IR + ) n_1 to the last (n — 1) components of Y(-) , using the 
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reflection matrix 



(2.1) 



1Z = I n _i — Q , where Q :- 





3 2 








\ 


4 























\o 





3n-l 





/ 



and I n _i is the unit (n — 1) x (n — 1) matrix. 

The main observation here, is that the spectral radius of the matrix Q is strictly 
less than 1. Indeed, the transpose Q T is an irreducible substochastic matrix, and can 
be made into a stochastic matrix by adding an absorbing point. Hence, by virtue of 
the Perron-Frobenius Theorem, the spectral radius of Q T , and hence also of Q , 
is strictly less than 1. 

All in all, we see that Theorem 1 of [13] is applicable here. We can now complete 
the definition of the process R(-) = (Rx(-), • • - , R n (-)) by imposing 



y^cfcj^o = c k R k 



k=l 



k=l 



(fll(-) - R 2 {-), • • • , Rn-l(-) ~ Rn(-)) = V HR (Rl(-) ~ M'), " " " > R n~l(-) ~ M')) ■ 

The strong uniqueness of the so-contructed process R(-) follows from the uniqueness 
of the solution to the mult i- dimensional Skorokhod reflection problem in Harri- 
son & Reiman (1981), Reiman (1984); see Theorem 1 of [T3]. 



(2.2 



Let us close this subsection by observing from (1.1) that the process of spacings 
Z(-) := (R 1 (-)-R 2 ( 



Ri 



R, 



i Rn- 



Rn{-)) 



for the process R(-) we just constructed, is a reflected Brownian motion (RBM) in the 



nonnegative orthant 
matrix 



\n-l 



with drift vector (b\ — b 



2, • • • , n -\ 



covariance 



(To 



(2.3) 



A 



-an 




\0 



- a 2 







-On 








and reflection matrix 1Z = I n _j — Q given as in (2.1). 



2.2. A boundary property of reflected Brownian motion. Throughout the 
paper, we rely many times on a generalization of the boundary property of reflected 
Brownian motion established in [22], which is of interest in its own right. 



Consider a continuous semimartingale Q(-) taking values in the orthant 



\n-l 



and satisfying 
(2.4) 



n-l 



(*)/ 



k=l 



where B(-) is an (n — l)-dimensional Brownian motion with a constant drift vector 
and a constant, nondegenerate diffusion matrix. Here, for each k £ {1, 2, . . . , n — 1} 
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and with 

'n — 1 



711 



k 



is an (n — 1) x m matrix with real entries; whereas 2)^(-) is a continuous 
(M + ) m — valued process, whose components are indexed by the sets Jc{l,2,...,n- 
1} with k elements, start at 2)^(0) = , are non-decreasing and satisfy 



POO 

(2.5) / Vl {Qj(s)>0} d2)J c) ( S ) = 0, JC {1,2,..., \J\=k, 
V < s < t : 2)f } (t) - Z) { l J 2 2l \s) < 2)f D (t) - ^ ll \s), 



Lemma 1. Let Q(-) be a process as in (2.4), and suppose that the matrix 



jeJ 

(2.6) 

for any J x C J2 C {1, 2, . . . , n — 1}. 

As we show in the following lemma, if the matrix is completely— S (see [22] for 
the definition and a characterization of completely— S matrices), then Q(-) can be 
identified with a reflected Brownian motion in the sense of HU, 1261 . 

is 

completely— S . Then all processes 2)^(-) , k — 2, 3, . . . , n — 1 are identically equal to 
zero. In particular, Q(-) is a reflected Brownian motion in the orthant (M + ) ?1 ~ 1 with 
reflection matrix 

Proof: The proof is similar to the proof of Theorem 1 in [22]. In view of Girsanov's 
Theorem, it suffices to consider the case that the drift vector of B(-) is equal to zero. 
We consider the functions <p E (•) , e £ (0, 1) as in the proof of Lemma 4 in [22]; these 
functions are harmonic for the generator of B(-) and bounded on compact subsets 
of (R + ) n_1 , uniformly over e 6 (0, 1) . We claim that, for any k = 1,2, ... ,n and 
any column v of the matrix D\( k \ there is a constant Ck < 00 depending only on the 
matrices 9^, 9^ and the diffusion matrix of B(-), such that 

(2.7) \/e > 0, x G (M+)"- 1 : v ■ V(f) £ (x) > -C h . 

Indeed, one can argue exactly as on page 94 of [22] in the derivation of the bound 
(24) there. Now, define the stopping times 

n-1 

r m = min(m,inf jt>0: \\Q(t) \\ + ^ ||2)W (t) \\ > m}) , m G N, 

fc=i 

where we wrote || • || for the Euclidean norm. Applying Ito's formula to the semi- 
martingale Q(-) and recalling from [22] that the functions £ (-) are harmonic with 
respect to the generator of B(-), one obtains for all e G (0, 1) and m G N: 



MQ(r m ))-MQ(0))= / V4> £ (Q(s))-dB(s) + J2Yl / vj-VMQ^dyf 

J ° fc=l Jea, - 70 



Here stands for the set of all subsets of {1,2,.. . ,n — 1} with elements and vj 
denotes the J-th column of 9 C {(' J '). Finally, taking the expectation on both sides and 
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using the bound (18) in [22] and the bound (2.7) above, one ends up with 
E[<f> £ (Q(r m ))) 



n-1 



: (g(0))]>-(loge + l)^ Ci E 



n-1 



k=l 



1 are defined as in 1 22 J and the 



where the positive constants Cj, fy, j = 1, 2, . . . , n 
constants k G {1, 2, . . . , n — 1} are as in (2.7). Dividing both sides of the latter 
inequality by (logs + 1) and taking the limit e I gives 



n-1 



lim > c.j E 



J'=l 







L{||Q( S )||< £ ^} dajj- 



(!)/ 



< 0. 



Thus, by Fatou's Lemma and the nonnegativity of the integrand, one can conclude 

,, j = 1,2,..., n-1, 



/ i {Q(s)=0} d2); 1) ( S ) = o, 

with probability one, and therefore also 
(2.8) 



poo 

/ l ww=o} dg)f(a) = 0, 



i = 1,2, 



n—1, 



with probability one. Finally, the backward induction argument of Lemma 5 in 
allows us to strengthen this statement to 



(2.9) 



poo 

Vj GJc{l,2,...,n-l}: / l {Qj{s)=0>jEJ} dZ)^ 

Jo 

almost surely. In view of (2.5) and ( |2.6 ), this finishes the proof. 



0. 



□ 



2.3. Absence of triple collisions. This subsection, and the one that follows, are 
devoted to the construction of the solution to the system of stochastic equations 
(1.4), when one of the following two conditions holds: 

[l-Qk)Ok > Ik a l+i » (! - Qk) °l ^ <?fc o-Li 



A 



k 



2,3, 



n — 1 



(B) q k 



1 + 



a 



fe-i 



fc+i 



2^ 



k = 2,3, 



n — 1. 



As we show below, each of these two conditions prevents collisions of three or more 
particles. It is not hard to see that neither of these two conditions implies the other. 

When all the collision parameters are equal to 1/2, as in (1.13), condition (B) 
mandates that the graph of the variances-by-rank k h- >■ a\ be linear; thus, under 
condition (B), Proposition [3] below generalizes the results in [18] on the absence 
of triple collisions to situations where particles feel local-time-like drag from their 
immediate neighbors, when they collide with each other. 

We start with a result ruling out triple collisions in the case n = 3 under the 
condition (2.10) below; in this three-dimensional case, condition (2.10) is weaker 
than each of the conditions (A) and (B). 
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Proposition 2. Suppose that n = 3, R±(0) — i?3(0) > and the condition 

(2.10) 2 <j\ > q 2 (a 2 2 + oj ) + q+ [a\ + a 2 2 ) 
holds. Then we have 

(2.11) P(3t>0: Ri(t) =R 2 (t) =Rs(t)) = 0; 
moreover, the converse is also true. 

Proof: We consider the reflected Brownian motion (-Ri(-) — R 2 {-), -R2O) — Rsi')) m 
the nonnegative quadrant, with the reflection matrix TZ = I 2 — Q of (2.1) and the 
covariance matrix A of (2.3) with n = 3 . Moreover, let O be the 2x2 orthogonal 
matrix such that £ := O' AO is a diagonal matrix, where the superscript ' denotes 
transposition. Then the process 

(2.12) Cr 1 O (R x {-) - R 2 (-), R 2 (-) - R 3 (-))' 

is a reflected Brownian motion in a wedge, in the sense of Varadhan & Williams 
(1984) in . Letting 

(2.13) V := diag(^) 

be the diagonal matrix whose diagonal entries coincide with those of A , and following 
the computations in subsection 3.2.1 of [18], one concludes that the normal vectors 
to the two sides of the wedge are given by the columns of the matrix 

(2.14) Af := C 1/2 OV- 1/2 ; 

whereas the reflection matrix of the new reflected Brownian motion in the sense of 
[24J is given by 

q -.= c- l,2 onv l i 2 -JT. 

Furthermore, as observed in the proof of Lemma 3.2 in [IB], the corner of the wedge 
is attainable by the new reflected Brownian motion, if and only if the sum of the two 
off-diagonal entries of the matrix M Q is nonnegative. 
Next, we note that 

(2.15) Jf'~Q = v-^nv 1 ' 2 -Jf'jT. 

Moreover, both off-diagonal entries of 7\F J7 are given by the negative cosine of the 
angle between the two sides of the wedge, which can be computed to be 

{A- X t 1 )'t 2 



((^- 1 e 1 )'c 1 ) 1 /2((^-i e2 )'e 2 ) 1 / 



where Ci, t 2 is the canonical basis of 1R 2 . Putting everything together, we can 
compute the sum of the two off-diagonal entries of the matrix AT Q as 



-92 / o ~ <?2 



'o\ + a 2 ^g 2 + a\ VK + ^X^ + ^I' 

Simplifying this expression, we obtain the condition ( |2.10 ) of the proposition. □ 
Now, we turn our attention to the case of general n > 3 . 
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Proposition 3. Suppose that n>3, and that condition (A) or condition (B) holds. 
Then no triple collisions are possible, that is, 

(2.16) P (3t > 0, l<i<j <k<n: R^t) = Rj{t) = R k (t)) = . 

The proof relies on an inductive argument and the following lemma. 
Lemma 4. Let n>3. Then, condition (A) is equivalent to the following condition: 

(2.17) Vl<j, k<n-l, j^k: [^(in-i - Q)] jk > 0. 

Proof of Lemma [4| We start by recalling that the entries of A -1 can be computed 
by the formula 



i-i _ r_i\j+k 



where A?' is the (n — 2) x (n — 2) submatrix of the symmetric matrix A , that one 
obtains by removing the j-th row and the k-th column from A. Next, we introduce 
the notation 



(2.19) ng := (n*2) (E^)> l<^i<^2<n, 

fc=fci fc=fci ™ 

and claim that the determinants det(.A J '' ) , l<j<k<n — 1 are given by 
(2.20) 

(-i)^det(^) = ^•••^(£-L)( £ JL) = ni( J] <K +1 . 

This claim can be verified easily, using induction over n and distinguishing the cases 
k > j + 2, k = j + 1 and = j. Therefore, for any fixed l<j<k<n— 1, the 
inequality 

(2.21) [^ 1 (ln-i-Q)] ifc >0 



of condition (2.17) can be rewritten in the equivalent form 

k fe-i fc+i 

(2.22) Hi ( af) IT£ +1 > q k U{ ( of) II£ + (1 - g fc ") ITj ( af) IT£ +2 . 

£=j+i £=i+i £=i+i 



In view of the strict positivity of the variances, the inequality (2.22) simplifies to 
(2.23) a\ W k+l > q- k U n k + (1 - q k ) a\ a 2 k+1 W k+2 . 

Now, we note the following relations among U k , H k+1 and H k+2 '■ 

(2-24) U k+l = U k+2 al +1 + a 2 k+2 ct| +3 ■ ■ ■ o%, 

(2-25) U k = U k+2 a\ a 2 k+l + (a 2 k+l + a 2 k ) cr 2 k+2 a 2 k+3 ■■■ a 2 n . 



Plugging these into (|2.23|) and simplifying further, we end up with 

T t+i ■ 



(2-26) > %°l 
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An analogous computation, now with 1 < k < j < n — 1, shows that the inequality 
(2.21) in this case is equivalent to 

(2-27) 0-qf)°i > 9t°J-i- 



All in all, we conclude that condition (2.17) is equivalent to condition (A). □ 

Proof of Proposition [3f We note first that, by the same application of Girsanov's 
Theorem as in subsection 2.2 of [18], or as in the proof of Lemmata 6 and 7 in [19], 



we need only consider the case b\ = b 2 = . . . = b n = in ( 1 . 1 ) . 

• We start by assuming that condition (A) is satisfied, and proceed by induction 
over n > 3 . For n — 3, very simple computation shows that condition (A) implies 
condition (2.10). The statement of the proposition for n = 3 follows then directly 
from Proposition [2j 

Next, we assume that n > 4 , and that the statement of the proposition holds under 
condition (A) for all v — 3, 4, . . . ,n — 1 (the "induction hypothesis"). Introducing 
the stopping times 

(2.28) = inf {t > : _ max (R k (t) - R k+l (t)) < 8 } , <5 > , 

k^ = 1,2. , ,7i 1 

we claim that 

(2.29) P (3 < t < t s , l<i<j<k<n: R^t) = Rj{t) = R k (t)) =0, 5 > 0. 
Indeed, for any fixed 5 > 0, the time interval [0,t$) can be written as 

M 

(2-30) [0,r,)= |J [rr,rr +1 ) 

m=0 

for some MgNU {oo} and with stopping times = r° < t} < . . . satisfying 

(2.31) 3k = k(m)e{l,2,...,n-l}: R k (t) - R k+1 (t) > 5 , t G [r 5 m , r™ +1 ) . 

Moreover, on each of the time intervals [t™,t™ +1 ) , < m < M, the system ( |1.1 ) 
splits into two subsystems of the same type which evolve independently, conditional 
on R^Tg 1 ), R 2 (rg l ), . . . , i^r™). Hence, (2.29) is a consequence of the induction 
hypothesis. 

In view of (2.29), in order to show (2.16) and complete the induction argument, it 
suffices to show 

(2.32) P (hmr 5 = oo) =1. 

To this end, we introduce the functions F , G : (!R + ) n_1 \{0} — > R + given by 

(2.33) F(z) := G(^) (3 - n)/2 , G(z) := {A^z.z) . 

Next, we apply the change of variable formula of Theorem 2 in [T3] to the function F 



and the reflected Brownian motion Z(-) formed by the spacings (see (2.2) and (1.9)), 
noting that the function F is harmonic with respect to the generator of the Brownian 
motion driving these spacings. We obtain the semimartingale decomposition 

F(Z(.)) = F(Z(0)) + M(.) + V(-), 
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where M(- 
and 



is a real- valued local martingale of the form M(-) = Y^l=i Jo £k(t) d/3 k (t) , 
V(-) =E f\d H -4d^ 1 -^d H , 1 )F{Z(s))dA.^ 1 \t) 



3 — n 



(G(Z(t))) 



(l-n)/2 



n-1 



vr 1 

A; (£-1) 



C + iA-(Vi))dA (V+1) W 



a process of finite variation on compact intervals. We also note that the processes 
£fc( • A 75) are all uniformly bounded, so the stopped local martingale M( ■ A r$) is 
in fact a martingale. 

In view of Lemma |4j condition (A) ensures that, in this last expression for the 
finite variation process V(-) , the coefficients appearing in front of the local time 
processes are all non-positive, so we conclude 

(2.34) V t > 0, 5 > : E [F(Z(t A t s ))\ < E [F(Z(0))\ . 

Moreover, if the limit lim^o t~s were finite with positive probability, then for any given 
real number c > there would exist real numbers t > sufficiently large, and 5 > 



sufficiently small, such that the value of the left-hand side in (2.34) would exceed c 



but this would then contradict (2.34). We conclude that (2.32) holds, and thus the 



proposition is established under condition (A). 

• Condition (B) simply paraphrases the skew-symmetry condition for the (n — 1)- 
dimensional Brownian motion Z(-) with reflection on the faces of the nonnegative 
orthant in the sense of Harrison & Williams (1987) (cf. [H], [15], [26] as well 
as section |2.5 below), and the result of Proposition [3] in this case can be found in 
Theorem 1.1 (in) of 126 1. □ 



2.4. Names. We can now combine the results of the previous two sections with those 



in [9], to give a construction of a strong solution to (1.4) satisfying (1.7) and to show 
its pathwise uniqueness. 

Theorem 5. Suppose that condition (A) or condition (B) is satisfied. Then, the 



system (1.4) has a strong solution satisfying (1.7), and such a solution is pathwise 
unique. 

Proof: We start with the proof of strong existence, which proceeds by induction 
over n. As remarked in subsection 1.3| for 



n 



2, the system of equations (1.4) is 



a special case of the system (4.13)-(4.14) in [5]. Therefore, we may deduce strong 
existence for n = 2 from their Theorem 4.2. Moreover, it is shown in that paper 
(see (6.13)-(6.15) in [9]) that the distributions of the random variables Xi(t) —X 2 (t), 
t > have a density with respect to the Lebesgue measure on R for every t e (0, 00) , 
and therefore 

(2.35) E[Ceb({t > : X 1 {t) = X 2 (t)})] = 

by Fubini's theorem. Hence, Xl(-), X 2 (-) satisfy ( fl~7| ). 



We now consider n > 3, and assume that a strong solution to the system (1.4) 



satisfying (1.7) has already been constructed for all v — 2, 3, . . . ,n — 1 and all choices 
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of drift, dispersion and collision parameters obeying condition (A) or (B) (the "in- 
duction hypothesis"). We shall construct a strong solution (Xl(-), X 2 (-), . . . ,X n (-)) 



of the system (|1.4|) consecutively on the random time intervals 
(2.36) 

where rj^t p , m 



[V2-P ,vl-p), H-p ,vl-p) , ■ ■ ■ , [Vz-p , Vz-p ) > P e N 



0,1, 



M(p), p G N are stopping times such that 



lim V ™ P = inf \t > : max (R%(t) - i?f +1 (t)) < 2~ p ) = r 2 - P , p G N 

tM(p) I fc=l,...,n— 1 J 



mfM (p) 

VpeN,m = 



M(p), 3 fc = k(p,m) : it£ +1 (i) - R*{t) > 2^, t 6 [v^V?^ 



and M(p) G N U {00}, p G N. We have recalled here the notation of ( 2.28| ), and the 
fact that the process of ranks (Rf(-), -Rjf (')> ■ ■ ■ , Rn(')) solves the system of equations 
(recall the discussion in section [T] for more details). 

On each interval [77™ j, , rf^ ) , we define (Xi(-), X 2 (•),..., X n (-)) by letting the 
processes (Xx(-), X 2 (-), . . . , X fc ( p m )(-)) evolve as a strong solution of the k(p, Tri- 
dimensional system corresponding to the first k(p,m) equations in (1.4), started 



at the point (Xi(ti™ p ), . . . , X k ^ m) {rf^ p )) ; and by letting (X k ^ m)+1 (-), . . . , X n (-)) 
evolve as a strong solution of the (n — k(p, m))-dimensional system corresponding to 
the last (n — k(p,m)) equations in (1.4), started at (Xk( p ,m)+i(vT-p) ■ ■ ■ ■ ^nijfz-p))- 
Note that the strong solutions to the lower-dimensional systems exist by the induc- 
tion hypothesis. The resulting process (Xi(-), X 2 (-), . . . ,X n (-)) is a strong solution 



to the system of equations (1.4), up to the random time 
(2.37) lim lim rf£- v = limr 2 - P . 

pt°° m\M (p) p\oo 



We can conclude from Proposition 3 that the quantity in (2.37) must be equal to 



infinity with probability one. Thus, we have constructed a strong solution to (1.4) 



for all t G [0, 00). Finally, it is clear that this solution satisfies (1.7) by the induction 
hypothesis and its construction. 

To prove pathwise uniqueness, we argue again by induction over n. For n — 2, 
pathwise uniqueness is a consequence of Theorem 4.2 in [5]. Now, let n > 3 and 
assume that pathwise uniqueness holds for all v — 2, 3, . . . , n — 1 and all choices of 
drift, dispersion and collision parameters satisfying condition (A) or condition (B) 
(the new induction hypothesis) . Next, suppose that X(-) — (Xi(-), X 2 , . . . , X n (-)) is 
another strong solution of ( 1.4[ ) satisfying (1.7), defined on the same probability space 
as the strong solution X(-) = (Xi(-), X 2 (-), . . . ,X„(-)) constructed above. Consid- 



ering the intervals in (2.36) consecutively and employing the induction hypothesis, 
we conclude that X(-) = X(-) must hold up to the time given by (2.37). However, 



as we have seen above, the latter must be infinite with probability one, by virtue of 
Proposition [3] This yields the desired pathwise uniqueness. □ 

2.5. Skew-symmetry and invariant measures. From Ichiba et AL. (2011) (see 
equation (5.8) in [20]), we know that the reflected Brownian motion of spacings in 
( J2T2] ) is skew-symmetric in the sense of Harrison & Williams (1987) (cf. [T4"] . 
[15]). if the following condition is satisfied 

(2.38) 2(V-A) = QV + VQ 
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Figure 1. 



Figure 2. 



with the notation V = diag(^4) of (2.13 ). Plugging in these equations the expressions 



for the matrices A and Q from (|2.3|) and (|2.1|), respectively, we can simplify this 
condition to 



(2.39) 



Qk = Ik 



a 



k-l 



fc+1 



2oi 



k = 2,3, 



1, 



that is, exactly the condition (B) 



In view of our assumption (|1.2|), the condition (|2.39|) amounts to the requirement 
(2.40) 



2o\ 



a 



2+°l 

2a\ 



2 erf 



Figure 1 shows the variances a k , 

lo 

k i- 



2 k 



1,2, 



n 



and Figure 2 the slopes 



log(or 2 



fe+i 



0"T 



log A; — log(fc — 1) 

of the function log A; h- > log(cr| +1 — af) , for n = 100 and a nonlinear choice of initial 
parameters of, a\ y a 2 , namely a\ = 0.1, a\ = 0.11, a 2 = 0.121. One can see that 
these slopes can be made to deviate significantly from 1 even when one only slightly 
perturbes a linear specification (all the slopes would be equal to 1 in a skew-symmetric 
specification of a model with symmetric collisions: = 1/2 , k — 1, 2, . . . , n — 1). 



Under the skew-symmetry condition (2.40), one can compute the invariant distri 



bution of the spacings process (2.2) explicitly, using Lemma 3.6 in the dissertation 
of Ichib A (2009) [TZ] ; one ends up with a product of exponential distributions with 
parameter vector 

(2.41) 7 = 2[dmg(A)]- 1 n- 1 (b 2 -b 1 ,b 3 -b 2 ,...,b n -b n _ 1 )' . 

Note that, by virtue of ( |1.2[ ) and ( 2.39[ ), the matrix TZ takes the form 



(2.42) 



TZ 



Ln-l 



Q 



I 


1 


-q 








\ 




-Q 


1 -(1-g) 















-(1-g) 1 


-Q 












-q 


1 






\ 












/ 
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where 

2 i 2 

(2-43) q: = q - = [1+ 



2 a? 



'2 

Remark 1. An interesting special case is the specification 

(2.44) bi — b 2 — ■ ■ ■ — & n -i = , b n = q n for some q > 0, 

which might be called a q-Atlas model by analogy with the term Atlas model intro- 
duced by Fernholz (2002) in [8] and studied further by Banner, Fernholz & 
Karatzas (2005) in [2]. Figure 3 shows on a log-log plot the capital distribution 
curve 

e R k (t) 

(2.45) k i — > , fc = l,2,...,100 

for such a g-Atlas model with n = 100, = 1, a\ — 0.1, a\ = 0.11, (T3 = 0.121, when 
the spacings process (-Ri(-) — R2{'), ^2(0 — -^3(')? • • • 5 -^n-i(') — -^n(')) takes the mean 
value under its stationary distribution. By comparing with the plots of real-world 
capital disribution curves from U.S. equity market data of the Center of Research 
in Securities Prices (CRSP) at the University of Chicago (see Figure 5.1 on page 95 
of [8]), one sees that a q- Atlas model can capture the concave shape of the capital 
distribution curve, as well as its linear structure at the top. □ 




0.01 - 



0.001 - 



10- 4 - 



10- 3 - 



Figure 3. 



We also note that, since the mapping 
(2.46) (b 2 - h, b 3 -b 2 ,...,b n - bn-!) 1— ► 7 

is bijective, one can determine the drifts up to an additive constant by fitting the 
vector 7 to the observed capital distribution curves, such as the ones in Figure 5.1 
on page 95 of [8]. 
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3. Scaling limit of asymmetrically colliding random walks 
Consider again the processes corresponding to the system of stochastic equations 



(1.1). Our objective in this section is to show that such processes with asymmet- 
ric local time components arise as scaling limits of random walks with asymmetric 
interactions upon "collision" . We start with the following informal description. 

Consider an n-dimensional continuous time jump process T(-) = (T^-), . . . , r n (-)) 
on the wedge 

H n = {(*!, ..,4)6Z": Zl > z 2 > . . . > z n } . 

The process registers the positions on the integer lattice of n particles that always 
maintain their order. We will label the particles by the elements of the set [n] = 
{1,2,..., n} where 1 refers to the rightmost particle and n refers to the leftmost. 

The movement of particles can be described as follows. Consider nonnegative 
parameters {ak,b k ,9 k ,9 k , k G [n]}. If there are no other particles at the current 
site of particle k, then this particle moves to the right at rate a k and to the left 
at rate b k , independently of every other particle. On the other hand, suppose that 
particles k, k + 1, . . . , k + £ are currently at the same location (a phenomenon we call 
"collision"), and this is the maximum length of the tie in the sense that the (k — l)-st 
and the {k + 1 + l)-st particles are not at that site. Then, at rate 9j}a k , particle k 
jumps to the right; and, at rate 9^ +e b k +i , particle {k + £) jumps to the left, while 
the particles in between do not move. 

We wish to take a diffusion limit of such systems. Before we go on, let us give 
a few examples which display the variety of behavior that we can expect from such 
particle systems. 

Example 1. Suppose each particle evolves according to a Poisson process with 
parameter a until a collision occurs. Thus every a k = a , every b k = , and the 
values of 9%, k G [n] are irrelevant. 

Consider, first, 9j} = 2, k G [n] . This corresponds via relabeling to the case when 
the Poisson particles are moving independently and are allowed to pass each other. 
On the other hand, when each 9j} = 1, the higher particle blocks the movement of 
the lower one when they collide. The resulting particle system then evolves as the 
well-known TASEP process (except that in TASEP collisions occur when particles 
are at adjacent sites instead of the same site). 

Example 2. Suppose the particles move according to i.i.d. simple symmetric random 
walks until collision. That is, a k = b k = 1, k G [n\. Moreover, we take 0% = 9j} = 9, 
k G [n] . Thus, when a particle is involved in a collision, this simply changes the rate 
of its next jump. The case 9 = 2 yields the ordered system of i.i.d. random walks. 
When 9 ~ 0, the particles get "sticky", while when 9 is very large, the particles can 
be thought of as repelling one another. 

Example 3. Consider, as before, particles moving according to i.i.d. simple sym- 
metric random walks until collision. That is, a k = b k = 1, k e [n]. Fix an i G [n], 
say i = 2. Choose the 9%, 9^ parameters to be much larger than the rest of the 
collision parameters. Thus, whenever particle 1 or particle 3 collides with particle 2, 
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the latter will almost immediately move away. Hence, particle 2 will remain invisible 
to its neighbors. 

Now, consider the triple collisions of particles 1, 2 and 3. Note that, with Q\ and 
$2 being very large, such triple collisions will happen for about the same duration of 
time as collisions between particle 1 and particle 3 in the particle system obtained 
by removing particle 2. In other words, if one removes particle 2, this should not 
influence the behavior of the other particles in a significant way. 

Example 4. We let the particles evolve according to i.i.d. standard Poisson pro- 
cesses until collision, as in Example 1. However, now we set 9^ = 0, with all other 
collision parameters being positive. Then particle 2 and particle 3 freeze forever the 
first time they collide. Thereafter, particle 1 moves as a standard Poisson process in- 
dependently from the rest of the particle system, while particles 4,5, ... ,n eventually 
all coalesce at the site of particle 2. 

Physical heuristics: To develop a feel for these processes, let us discuss briefly the 
mechanics involved in these collisions. In classical models of collision (see [12]), the 
particles behave as hard billiard balls of infinitesimally small radius. This models 
elastic collision which we now explain. 

Suppose two particles collide on the line. Particle 1, on the right, has mass mi 
and velocity U\ right before the collision. Particle 2, coming from the left, has mass 
m 2 and velocity u 2 - in elastic collisions the momentum and the kinetic energy are 
preserved before and after the collision. These two preserved quantities uniquely 
determine the velocities V\ and v-i of the two particles right after the collision, i.e., 

U\{m\ — 777.2) + 27772^2 M2( m 2 — m l) + 2777i«i 

Vi = , v 2 = . 

777i + 1712 mi + n%2 

Consider again our jump processes. We follow the heuristics of [12J and assume 
that particles jump together at discrete times. A similar but slightly more lengthy 
analysis can be done for continuous time, when we let one particle stay still and get hit 
by the other particle. We consider the Einsteinian viewpoint, that the n particles are 
bombarded on all sides by other small particles which lead to their random motions. 
Suppose particles k and k + 1 are adjacent and are of masses mi and m 2 respectively. 
We will think of the rate of jumps as the speed in the appropriate direction. 

In the next small time interval, these two particles either do not collide, and move 
away from one another at velocity and — b k +i respectively; or they collide, and 
ui = —bk , U2 = cifc+i • Thus, if the collision is elastic, we will observe that the 
velocity (total jump rate to the right) for particle k is 

-6 fc (mi - m 2 ) + 2m 2 a k+ i 
Ok '■— ctk H ; , 

777! + 777 2 

and the total jump rate to the left for particle (k + 1) is 

a k+ i{m 2 - mi) - 2mib k 

Ok+i '■— —Ok+i H ; ■ 

m i + m 2 

Clearly, unless we impose specific constraints on the parameters, these are not equal 
to 9j}ak and —6^ +1 b k+ i , respectively. Hence, these collisions are not elastic in gen- 
eral. 
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Certain special cases are worth mentioning. Suppose all masses are equal to one. 
Then, under elastic collision, 5k = + ctk+i and 5k+i = —bk — 6fe+i. This is the case 
when particles exchange velocities and can be thought of, via relabeling, as crossing 
over. A specific choice of 9j} and 9% +1 would capture this scenario as in the case of 
9% = 2 in Example 1 and 9 = 2 in Example 2 above. 

On the other hand, suppose that all a>k and bk are equal to 1. Then, under elastic 
collision, 9j} = 5k = 4m2/ (mi + 777,2) and 9% +1 = 5k+i = 4mi/(mi+m 2 ). In particular, 
if one of 9j} + 0% +1 is greater (or, less) than 4, the mechanics generates excess energy 
(or absorbs energy) that cannot be explained by elastic collision. 



Under suitable assumptions on the parameters, we establish in Theorem 11 diffu- 



sion limits similar to (1.1) for these particle systems. There is an apparent paradox 
here. Consider that diffusion limit for a = b = aj^ = a^ = l for all k ; somewhat 
surprisingly, the limit turns out to depend only on the ratio 9f} l®k+\ ■ In other words, 
the diffusion limit is the same whether Of} + 9% +1 is equal to 4 or not. The collisions 
among the limiting diffusion particles are always elastic, and the quantities g^T's and 
q^s can be thought of as the proportions of total mass shared by the colliding par- 
ticles. This is a consequence of the fact that the occupation time of collisions has 
Lebesgue measure zero in the limit, as will be made clear in the proof. 

To get a true inelastic limit, one has to let 9^, 9^ go to zero suitably in the diffusion 
scaling; then, presumably, one would obtain sticky colliding Brownian particles. 

3.1. The modified Skorokhod problem. We start with a set of parameters: 
a,b, (X^,a^,9^, k G [77]), (Af , af* , 9j} , k 6 [77]). At this point we only assume that 
(T^ , aj} are stricly positive and 9%, 9j} are nonnegative for every k. 

Taking a diffusion limit requires considering a sequence of interacting jump pro- 
cesses as described above. We shall generalize the setup by allowing non-exponential 
waiting times for the jumps. Let the sequence of interacting jump processes be in- 
dexed by iV > 0. We fix a value of iV and a probability space rich enough to support 
mutually independent sequences of i.i.d. random variables (u^(i), i G N), k 6 [77] 
and (u^ii), i G N), k 6 [n], all taking only positive values. These random variables 
denote the inter-jump times of the particles (the superscripts "L" and "i?" standing 
for leftward and rightward jumps, respectively). Assume that, for any fixed k G [n], 



E H(l)] = ( b + A=) , VarK(l)] = (^) 2 , 

a) v v -y , 

E K(l)] = (a + ^L) , VarK(l)] = K) 2 - 

Next, we define the corresponding partial sum processes 

3 

U k L (0)=0, ut(j) = X>£(0, JGN, ke[n], 

i=l 

f/ fc R (0)=0, l/ftf) = Z>*W' J6N, ke[n], 



i=l 



20 



IOANNIS KARATZAS, SOUMIK PAL, AND MYKHAYLO SHKOLNIKOV 



and the corresponding renewal processes 

S^t) = max {j > : U^j) < t} , t>0, k E [n], 
S*(t) = max {j > : U*(j) < t} , t > 0, fc G [n]. 

Finally, we denote by (J-" t , t > 0) the filtration generated by the processes S k (-), 
S k (-), k E [n). Informally, for each k E [n], the process S'l'(-) (resp., S k (-)) records 
the leftward (respectively, rightward) jumps of the A;-th particle from the right, as 
long as this particle is not involved in a collision. 

To describe the effect of collisions we shall use a stochastic time change. The 
following lemma encapsulates the idea that the leftward (respectively, rightward) 
movement for the of particle k can either be blocked, or proceed at a different rate, 
depending on whether there is a collision with particle {k + 1) (resp., particle (&; — 1)). 
Such results are standard in the Queueing Theory literature (e.g., section 2 in the 
seminal article [21] by Reiman (1984)), so we omit the proof. 

Lemma 6. For every N E N, and any (71,72, • • • ,7«) G M n , there exists a system 
of jump processes T(-) = (r fc (-), k G [n]) taking values in M n and progressively 
measurable with respect to T t , t > 0, that satisfies the following set of equations 
pathwise: 

T k (t) = lk - S L k {T k L (t)) + S* {T k R (t)) , where 

(3.2) TfcL ^ = J 1 {Qk-i(s)>o] dT fe (s) + d k 1{0*_iW=o} dT fc (s), and 

T k( t )= 1 {Q fc (s)>o} dr fc _i(s) + 6 k / l{Q fc(s ) =0 } dT fe _i(s), 
Jo Jo 

k G [n] . We have denoted here by 

(3.3) Q(.) = (Q fc (.), ke[n-l])= {T k (-) - T k+1 (-), k G [n - 1]) 
t/ie process of gaps, and have set 

T k (t)= l {Qfe(s)>0 }ds, fc = 0,l,...,n, 
^0 

mi/i t/ie convention 1{q o( . )>0 } = 1{q„(.)>o} = 1- 

It should be noted that, although we have suppressed the index iV from the no- 
tation in the lemma above, this parameter determines the drifts of the coordinate 
processes. The key to passing to the scaling limit is to understand the time-changes 
involved. Our strategy is the following: (i) express the process of gaps as a Skorokhod 
map applied to a suitable "noise process" ; (ii) show that the sequence of distributions 
of gaps is tight; (iii) and finally, show that tightness implies the convergence of the 
appropriately rescaled process T(-) to a semimartingale of the type described in (1.1 ). 

We start by analyzing the process of gaps. To this end, we define the centered 
processes 

S k (t) = ^ L (t) - bt, t > 0, S k (t) = S*(t) -at,t>0 



(3.4) 
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for all k G [n\. Also, we define the following processes measuring the time spent in 
the various collisions: 

h(t)=t-T k (t), t>0, fc = 0,l,...,n, 

h,k+i(t) = / l{Q fc ( s )=Q fc+1 ( s )=o}ds, k = 0, 1, . . . ,n - 1. 
Jo 

With this notation, 

Q*(f) = 7. - 7Hi - ^ (T fe L (t)) + S« (Tf(t)) + S[ +1 (T t V(*)) - (T^(«)) 
= Ik- 7* + i - Si (T k L (t)) + S« (T R (t)) + S L k+l (T£ +l (t)) - S« +1 (T« +1 (t)) 
-b {T k L (t) - T k L +1 (t)) + a (T k R (t) - T k R +1 (t)) , k G [n - 1]. 

Moreover, 

Tkif) = T kif) + {Ok ~ 1) / 1 {Q k -i(s)=o,Q k (s)>o} ds 

J 

= i - 7 fc (i) + - 1) (i^t) - jf l Wfc _ l(s)=Qfc(s)= o } ds) 

= t - j fc (t) + - f) - (9 k - 1) h- ltk (t), t > o 

and 

(3.5) T R (t) = t - h^{t) + (9 R - l) I k {t) - (9 R - l) 7 fc _i, fc (0, * > 
for all k G [n\. 

Next, we introduce the process X(-) = G [n — 1]) via 

**(•) = -S L k {T k L (.))+S* (T*(.)) (T^ +1 (-)) -Ci Pfti(O) , 

G [n]; this will play the role of "noise process", to which a Skorokhod map will be 
applied to obtain the process of gaps Q(-) = (<5fc(-)> ^ e [ n ~~ !])• 

Combining everything so far, we obtain the representation for these gaps 

Qk (t) = Ik ~ lk+i + X k (t) 

- b (I k+1 (t) - 4(0) - b (0 k - 1) I k ^{t) + b (0 k - 1) J fc _ 1>fc (f) 
(3-6) + b (9 L k+1 - 1) I k {t) - b (9 L k+1 - 1) I kM1 (t) 

+ a (I k (t) - 7 fe _i(0) + a (9 R - 1) J fc (t) - a (9 R - l) I k ^ k (t) 

- a (9 R +1 - 1) I k+1 (t) + a - 1) I k , k+1 (t), k G [n]. 

Now, it is very convenient to write this equation in a more transparent, matrix- 
vector notation, so we introduce suitable "reflection matrices" . Let £H be the (n — 
1) x {n — 1) matrix such that, for every k G [n — 1], all entries in the k-th row of 
are zero, except the (k — l)-st, the k-th and the (k + l)-st. These are given by 

( o 7) „ _ a + „ » _ & + a(g fc fl + i-i) 
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Note that the columns of the matrix I n _i — y\ add up to one. Similarly, define the 
(n — 1) x {n— 2) matrix %\ such that, for each k G [n — 1], all entries in the k-ih row 
of are zero, except for the (k — l)-st and the fc-th, which are given as 

(3.8) = -a (8* - 1) + b (9 L k - 1) , &m = a (9* +1 - l) - b (6 L k+1 - l) . 
Finally, we introduce the stochastic processes 

(3.9) 2) fc (.) = (a0* + 60£ +1 ) fe G [n - 1], ?j*(0 = W(-), fc e [n - 2]. 



With this notation, we have the following cleaner matrix- vector analogue of (3.6): 

(3.10) Q(-) = Q(0) + X(-) + <K2)(-) + 

We shall refer to this equation as the modified Skorokhod representation. 

3.2. The diffusion limit. To be able to pass to the diffusion limit, we make the 
following assumptions on the parameters. 

Assumption 1. Assume that a > 0, b > 0, 
(i) the entries of 91 satisfy 

a + b(9 L k -l) rn & + a (0* i) 



L 

fe+2 



(m^ i/iere is a fco £ + 1] suc/i t/iat £/&e numbers a (6 k — l) —b (6 k — lV k G [n] 

are nonpositive for all k < k and nonnegative for all k > k ; 
(Hi) for some e > 0, we have 

sup max (E [4(l) 2+£ ] + E [u* (1) 2+£ 1 ) < oo. 
jv e N fce[n] 

A simple case, in which parts (i) and (ii) of Assumption [I] are satisfied, is given 
by ^fe — @k > 1, k G [n]. Next, for each m 6 N, we write D m [0, oo) for the space of 
right-continuous IR m -valued functions on [0, oo) having left limits, endowed with the 
topology of uniform convergence on compact sets. 

On the strength of part (iii) of the above assumption, the following result follows 
from Theorem 14.6 in [5] (note that one can improve the topology in the conclusion 
of Lemma [7] to the locally uniform topology, by noting the path continuity of the 
limit process, modifying the paths of the jump processes to continuous, piecewise 
linear functions, and using the fact that the Skorokhod topology relativized to the 
space of continuous functions coincides with the locally uniform topology there). 

Lemma 7. The distribution of the process 

(J=S L k (Nt),t>0, J=sf(J\rt),t>0, fce[n]) 



converges weakly in D 2n [0,oo) to the law of the vector of independent processes 
(Z^(-), Z^(-), k G [n]) . Here, for each k G [n], Z k {-) is a Brownian motion with 
drift coefficient \ k and diffusion coefficient b 3 ^ 2 a k , while Z k (-) is a Brownian mo- 
tion with drift coefficient X k and diffusion coefficient 
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Next, we introduce the rescaled versions of the gap process by 

1 
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Q N (-) = k€[n-l]) 



and define 2) (•), 2J (•) accordingly. Then, 



Qk (Nt) , t > 0, fc e [n - 1] 



(3.11) 

where 
(3.12) 

x N k {t) 



Q 



N i 



Q N (0) +X (■) + 2}2T(-) + 9*2J 



a/TV L 



^ (T fc L (iVt)) - S n k {T k R (Nt)) ~ S; +1 (T k L +1 (Nt)) + Sti 



t > 0, k G [n — 1]. We can state now our main limit theorem for the process of gaps. 



Theorem 8. Suppose Assumption [I] holds, and that \im^^ oo Q N (0) = £(0) m cfe- 
tribution for some random vector £(0). Lei £(•) fre a reflected Brownian motion in 
(M + ) n_1 initial condition £(0), dn/t vector 



(3.13) 



&=(-A£ + A* + A£ +1 -A* +1 , fce[n-l]) 



diffusion matrix 21 = (2tfc^)i<fc^<n-i given fry 
(3.14) 

H 3 K + i) 2 -^ 3 K L +1 ) 2 



2t*,i 



if 
if 



fc + l, fc e [n-2], 
k-1, £e [n-2], 



a*(aZy + a*(<T« +i y + bZ(o-£y + b*(a% +i y if k = le[n-l], 



otherwise, 



and reflection matrix 21 as in ( |3.7 ). 

Then, the processes Q N {-), N > converge in distribution to £(•), as iV — >■ oo, in 
.D n-1 [0,oo). Moreover, the processes 



(?)?(•). * e [n - 1]) = « + 6^ L ) / 1 {0 * W=0) da, fc e [n - 1] 



(3.15) 



converge in distribution in D n 1 [0, oo) to the process of local times accumulated by 
£(•) on the respective faces of the orthant (IR + ) n ~ 1 ; and the processes 



(3.16) 



2jf(-), fee [n-2] 



L WfM=Qf +1 (s)=o} 



ds, k E [n - 2] 



tend to zero m distribution in D n [0, oo 



We remark at this point that the limit process can be viewed as the process of gaps 
for a semimartingale as in upon the appropriate identification of parameters 
(see Theorem 11 below for the details). The proof of Theorem M relies heavily on 



an extension of the invariance principle of Williams [27], which we introduce and 
prove in the following subsection. 
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3.3. Oscillation estimates for modified Skorokhod maps. In this subsection, 
we consider families (Q N (-),X (•), 2) iV (-), 2) 7V (-)), iV > of processes with right- 
continuous paths having left limits, which satisfy the following properties. The pro- 
cesses Q N (-), N > take values in the orthant (IR + ) n_1 ; the components of the 
processes %) N (-), N > and %) N (-), N > are non-decreasing and start at 0; the 
modified Skorokhod representation 

(3.17) q n (-) = x N (-) + m<r> N (-) + mz) N (-) 

holds for each N > 0, where 9^ and *R are as in ( |3.7 ) and ( |3.8 ), respectively; and the 
following are true: 

POO 

(3.18) / l {Qns)>0} d2)f( S ) = 0, ke[n-l], 
Jo 

/oo 
(l Ws>0} + l W & lW >o}) d?)?(s) =0, k G [n - 2], 

(3.20) V0< S <t: ^(t)-^( S )<min(^(t)-^( S ),^ +1 (t)-^ +1 ( S )). 

In this situation, we have the following extension of the invariance principle of [27J. 

Proposition 9. Under parts (i) and (ii) of Assumption^ suppose that the initial 
conditions Q N (0) = X N (0), N > converge in distribution to a random vector £(0). 
Then: 

(a) If the family of processes X N '(•), N > is tight on Z^ n_1 [0, oo), then the same is 
true for (Q N {-), X~ N (-), Z) N (-), ?)*(■)) , N > on D An ~ 5 [0, oo). 

(b) In the situation of part (a), any limit point (Q°°(-), X°°(-), 2)°°(-), ^t/j 
continuous paths satisfies 

(i) Q°°(t) = X°°(t) + <RZ)°°(t) + <Kjj°°(t) G (M+) n " 1 ; t > 0. 
(ii) The components of 2}°°(-) and 2)°°(-) are non-decreasing, and satisfy 2)°°(0) 

0, 2j°°(0) = wift probability 1. 
fmj TTie following two identities hold with probability 1 : 

/•oo 

(3.21) / l {Qr(s)>o} d2)-00 = O, fce[n-l], 

JO 

/>oo 

(3-22) / (l {Qr(s)>0} + l {or (s)>0} )d2)^( S ) = 0, fce[n-2]. 

Jo 

The proof will follow the ideas in [27], with additional complications caused by the 



last summand on the right-hand side of (3.17). As there, for any m G N and any 
function / G D m [0, oo), we introduce the notation 

(3.23) OscgCf) := sup max \f k (t) - f k (s)\. 

The proof of Proposition [9] is based on the following lemma. 

Lemma 10. Suppose that the functions q(-), sc(-), t)(-) G D" _1 [0, oo), t)(-) G D n_2 [0, oo) 
the analogue 
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of (3.17). Moreover, suppose that the function q(-) takes values in (I 



m-l 



, 0(0) = o, 



t)(0) = ; all components of [)(■), rj(-) are non- decreasing and 

POO 

(3.24) / l {gfc(s)>0} dt) fc (s) =0, ke[n-l], 
Jo 

POO 

(3.25) / (l {gfc(s )>o} + l{ gk+1 (s)>o}) dQfc(s) = 0, k E [n — 2] , 
Jo 

(3.26) V0 < s < t : tj fc (t) - < min(tj fc (t) - t) h (s), t) k+l (t) - t) fe +i(s)). 

Then, there is a constant C G (0, oo) depending only on the entries of 91 and ¥i (but 
not on the particular functions q(-), x(-), t)(-), t)(-)) such that 



(3.27) 



V < h < t 2 : Oscg (g) + Oscg (t)) + Osc*J (t)) < C Osc£ (x 



Proof of Lemma 10: We proceed by induction over the dimension n > 2. For 
n = 2, the lemma is a direct consequence of Theorem 5.1 in [27], since tj(-) = in 
this case. From now on, we take n > 3 and assume that the lemma holds for all 
v = 2, . . . , n — 1. Consider first the case that there exists a k e [n — 1] such that 
Ofc(^) = t)fc(*i). Then, (|3.26|) shows that tj k (t 2 ) = %(h) and tj k -i(t 2 ) = Qfc-i(ti). 



Therefore, the induction hypothesis and the same argument as on pages 15-16 in [27] 
imply together that (3.27) holds in this case. 



m— 1 



such that 



We now claim the existence of a vector A = (Ai, • • • , A n ) G (I 
(A'D^fc > 1 and (A'D^fc > for all k G [n — 1], where A' stands for the transpose of A. 
Indeed, noting that the off-diagonal elements in every column of 9t are negative and 
add up to —1 (see Assumption |l[i)), and that the entries of each column of add 
up to and satisfy AssumptionTjjn), we see that we may choose Ai, A 2 , . . . , A n _i as 
/(ai), f(a 2 ), • • • , /(a n -i) f° r appropriate numbers = a x < a 2 < . . . < a n ^i < 1 and 
an appropriate strictly concave function / : [0, 1] — > (0, oo). Next, we fix such a A 
and proceed as in the derivation of the inequality (18) in [27J to conclude 



'n-l 



(3.28) 



V0 < h < t 2 : Oscg(t)) < X'(q(t 2 ) - q(h)) +\J2 Xk ) 0sc ^^)- 



vfc = l 



Recalling (3.26) and proceeding as in the derivation of the estimate (20) on page 17 



of [27], we can find a constant C < oo such that 

(3.29) V0 < h < t 2 : Osc^(g) < cio^Mx) + max q k (t 2 ] 

V fcsfn— 11 



The lemma can be now obtained from (3.29) in the same manner as Theorem 5.1 



in [27] is obtained from the estimate (20) there. For the convenience of the reader, we 
give a short sketch of the proof. Fix < ti < t 2 and let < K < oo be a constant, 
whose value will be determined later. We distinguish between two cases: 



(a) 
(b) 



Qk{ti) > i^Osc^(x) for some k G [n — 1]. 
qh(h) < KOstffx) 



for all k G [n — 1] . 



In case (a), let r be the first time t G [ii,^] that q k (t) = and set r = oo if the 
latter event does not occur. If r = oo, then (3.27) holds by the argument in the first 
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paragraph of this proof as a consequence of the induction hypothesis. Now, suppose 
that r 7^ oo. Then, the same argument based on the induction hypothesis shows that 

(3.30) Osc^g) < COsc T tl (x), 

where we have increased the value of C < oo if necessary. Next, let C < K < oo. 
Since we are in case (a), we have 

(3.31) q k {r) > q k (h) - Osc£(g) > (K - C) Osc^x) 

which is a contradiction to r ^ oo. 

In case (b), we distinguish two possibilities: 

(i) q k {t) < KOsc\\(x) for all k G [n - 1] and t G [t u t 2 ]. 

(ii) Qkit) > KOsc t 2 (x) for some k G [n — 1] and t G [ti, t 2 ] • 



> 0, 



In case (i), the inequality (3.27) follows from 
(3.32) 



Osc^(g) < sup max q k {t), 
ti<t<t 2 fce[n-i] 



(3.28) and (3.26). In case (ii), we let r be the first time t G [£1,^2] such that q k (t) > 
KOsc^ix) for some k G [n — 1]. Then, we split the time interval [ti, ^2] into [ti> T ] 
and [r, t 2 ], and argue as in case (i) on [ti,r] and as in case (a) on [r,t 2 ]. □ 

Proof of Proposition |9| First, we claim that part (a) of Proposition [9] is a con- 
sequence of the inequality ( |3.27 ) in Lemma 10 Indeed, the necessary and sufficient 

conditions (a) and (b) of Corollary 3.7.4 in [6] hold for any subsequence of X (•), 
N > 0, and carry over to the same subsequence of (Q N (-),X (-),Z) N (■)), 

N > via the inequality (3.27). For more details on the same argument, please see 



the proof of Theorem 4.1 in 

Now, we turn to the proof of part (b) of Proposition [9 and let (Q°°(-), X°°(-), 2)°°(-), 
2}°°(-)) be a limit point as there. The properties (i) and (ii) for it can be seen by using 
the Skorokhod Representation Theorem in the form of Theorem 3.1.8 in [6] (noting 

the path continuity of the limit) for the subsequence of (Q N (-),X (•), 2) JV (-), 2) iV (-)), 
N > 0, which converges to that limit point, and taking the almost sure limit on both 
sides of the identity ( 3.17| ) through this particular subsequence. To deduce property 
(iii), we let g : [0, 00) — > [0, 1] be a continuous function such that, for some 5 > 0, it 
holds g(a) = whenever < a < 5, and g(a) = 1 whenever a > 28. By arguing as in 
the second half of the proof of Theorem 4.1 in [27], we conclude that the quantities 

POO 

(3.33) / g (Qf (s)) d2)f (s) = 0, ke[n-l], 
Jo 

POO 

(3.34) / g(Q»(s))+g(Q" +1 (s))dZ)»(s) = 0, k G [n - 2] 
converge to the quantities 

POO 

(3.35) / g{Q?( 8 ))dSQ?(8), k E [n - 1], 



(3.36) 



9 (Qk(s)) + g (Q&M) d2TO, k G [n - 2] 
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when we pass to the limit through the same subsequence as before. Now, letting 
5 l , we obtain the property (hi). □ 

We are now ready for the proof of Theorem [8] 

Proof of Theorem 8| Step 1. Consider the family of processes X N (■), N > 

in gig ). The family of processes ^S^Nt), t > 0, ^sf{Nt), t > 0, k G [n]j, 

N > 0, without the time change, is tight by Lemma [7} Moreover, the time-changes 
in (3.12) are Lipschitz functions of time. That is, there exists a constant O < oo 
such that 

V0 < a < t : max (T k L (t) - T k L (a) , T* (t) - T k R (a)) < 0(t -a), k G [n]. 

Hence, using the necessary and sufficient conditions for tightness of Corollary 3.7.4 
in [6], the tightness of X (•), N > is easily verified. However, we still have to 
identify the limit points. 

Step 2. At this stage, we can use Proposition [9] to conclude that the family 
^Q N (-),X N (-),ty N (-),Z) N (-)^, N > is tight. Now, recall that, for any N > 0, 

2) Ar ( - ) and 2} JV (-) can be expressed, as in (3.9), in terms of the times that the par- 
ticles spend in collisions. The tightness of %) N (-), N > and %) N (-), N > now 
shows that the processes 

(3.37) —I k (Nt),t>0, ke[n-l], —I ktk+1 (Nt),t>0, k e [n - 2] 
all tend to zero in D[0, oo). 



Step 3. Putting together the conclusion of Step 2, (3.4) and (3.5), we conclude that 
each of the processes 

±T£ (Nt) , t > 0, ke[n], ^T k R (Nt) , t > 0, k G [n] 

converge to the process t, t > in _D[0,oo). With the help of the Lemma on page 
151 in [5], preceding Theorem 14.4, we deduce f that the joint distributions of 

1 St (Tt(Nt)), t > 0, -±=s2(T*(Nt)) t t>0,ke[n] 



'N VN 

now with the time-change, converge on D 2n [0, oo) to the limiting distribution de- 
scribed in Lemma [7] (note that we can improve the topology used in [5] to the 
topology of uniform convergence on compacts by observing the path continuity of 
the limit process, modifying the paths of the jump processes to continuous, piecewise 
linear paths and using the fact that the Skorokhod topology relativized to the space 
of continuous functions coincides with the locally uniform topology there). 

Step 4. Step 3 implies that, in the limit A^ — >■ oo, the processes X N '(■), A^ > con- 
verge vaguely in D n_1 [0, oo) to a multidimensional Brownian motion with drift and 
diffusion coefficients as described in Theorem [HJ To conclude the proof, we note that 
every limit point X°°(-), 2)°°(-), S°°(0) of the collection (Q N (-),X N (-)^ N (■). 

2} JV (-)) , N > has continuous paths by (3.27) and the fact that the processes X (■), 
A^ > converge to a process with continuous paths. Thus, by part (b) of Proposition 
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|9| every such limit point satisfies the properties (i), (ii), (iii) there, with X°°(-) being 
the multidimensional Brownian motion just described. Lastly, by Lemma [TJ we can 
identify Q°°(-) with £(•), 2)°°(-) with the boundary local times of £(•), and deduce 
that 2)°°(-) = 0. This completes the proof. □ 

Finally, we consider the limit of the entire collection of jump processes (!"&(•), k G [n]) 
(we refer the reader to the expression (3.2)). Let T N (-) denote the vector of centered 



and scaled jump processes given, for every k G [n], by 

ri(t) := 4=r*(0) + -±= [S« {T k R (Nt)) - S L k {T k L (Nt))} - (a - b)tVN, t > . 
V iV v iv 



Theorem 11. Suppose that Assumption ll\ holds and that linijv->oo iV -1 / 2 rfc(0) = 
R(0) G W n in distribution. Further, let R(-) = (Ri(-), R 2 (-), . . . ,R n (-)) denote the 
continuous n- dimensional semimartingale taking values in W n and satisfying 

R k (t) = ^(0) + (Af-A9t+(a 3 (af) 2 + & 3 K) 2 ) 1/2 /3,W 



a 



k G [n] w'^A toe same notation as in (1.1). Then, the processes T (•), iV > converge 



in D n [0, oo) to toe process R(-) described above in the limit N — > oo. 
Proof. The main observation is that, for any fixed iV > 0, we have 



1 _ , . 1 



(3.38) i a ^~ 1)+ V m b ^- 1 ) +a ^N (t) 

+ (-a (0* - 1) + 6 (C{ - 1)) 2)f-i(0, t > 0, 
G [n] . The same steps as in the proof of Theorem [8] now show that the processes 



in the first line of (3.38) converge jointly, in distribution, to the components of 
a multidimensional Brownian motion with drift and diffusion coefficients as in the 
statement of this theorem, the process 2) iV (-) converges to the process of boundary 
local times of a reflected Brownian motion as in Theorem [8j whereas the process 
%) N (-) converges to zero. Therefore, the processes T Ar (-), N > must also converge 
in distribution, and one can identify the limit of the "noise part" with the appropriate 
multidimensional Brownian motion and the limit of the "local time part" with the 
local time part in the decomposition of the process R(-) in the statement of the 
theorem. This completes the argument. □ 

4. Additional Determinantal Structures 

This last section studies conditions on the parameters b±, b 2 , ■ ■ ■ , b n , o~i, a 2 , . . . , a n 
and qf, qf, . . . , q^ , under which the process R(-) = (Ri(-), R2{-), • • • > Rn( - )) of ranks 
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as in (1.1) has a probabilistic structure of determinantal type, in the sense that its 
transition densities are of the generalized Karlin-McGregor form 

n 

(4.1) p(t,r,f) = ** II f k ^ k \^^)-r k ) . 

cr& n k=l 

Here S n is the group of permutations of a set with n elements, whereas K a , <r G 
§„ and / fc ^, 1 < k, £ < n are suitable real numbers and real-valued functions, 
respectively. 

This question is motivated by the following two extreme cases: 

(i) If one considers b\ = b 2 = ■ ■ ■ = b n , 0\ = a 2 = . . . = a n and qf = q 2 = ■ ■ ■ = 
q^ = 1/2, then the ranks evolve as the ordered system of n Brownian motions, each 
with drift b\ and dispersion o\. Therefore, in this case 

n 

(4.2) P(t,r,r) = JJ ip bu(T1 (t, f a{k) - r k ) , 

cr£S n k=l 

where ipbi,oi(t, •) denotes the Gaussian density with mean b\t and variance aft. 

(ii) Now, consider the case of b\ — b 2 — ■ ■ ■ — b n — , a\ — a 2 — . . . — a n — 1 and 

(4.3) g+ = l, fc = l,2,...,n, = 0, fc = 1, 2, . . . , n . 

Then, the process of ranks is given by the continuous version of the totally asymmetric 
simple exclusion process (TASEP) treated in detail by Warren (2007) in section 4 
of [25]. In this case, the transition probability densities are of the form 

n 

(4.4) p(t, r, f) = ]T J] (t, fv(*) - r fc ) 

for suitable functions ?/' fe ' £ , 1 < k, £ < n (see Proposition 8 in [25]) and with sgn(cr) 
standing for the signum of a permutation cr 6 §„. 

The main result of this section is the following proposition. 

Proposition 12. Suppose that, for every e > 0, the transition probability densities 
of the process in (1.1) belong to the function space C;,((e, oo) x W n x W n ), and are 
continuously differentiable in the first coordinate and twice continuously differentiable 
in the second coordinate with derivatives in Cb((e, oo) x W n x W ra ). 



Then these transition densities are given by (4.1) with suitable real constants n a , 
cr G S n and real-valued functions fk,e , 1 < k, I < n , if and only if: b\ — b 2 — ■ ■ ■ — 
b n , o\ = o 2 = . . . = o n and for each k — 1, 2, . . . , n — 1 we have, either 

(i) Qk_=<lk+i = V 2 , or 

(ii) q^ = 0, = 1 , or 
(Hi) ql = I, qt +1 = 0. 



Moreover, if this is the case, one may choose K a = 1 , cr G S n and the functions f 
1 < k, £ < n according to the formulas (4.17)-(4.20) below. 
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Proof. We pick a bounded and continuous function g : W n — > R vanishing in a 
neighborhood of the boundary of W n , an arbitrary real number T > , and consider 
the martingale 



(4.5) 



/ p(T — t, R(t),r) g(f) df =: F(T — t, R(t)), 0<t<T. 



Since the change of variables formula of [T3] is applicable to the spacings process Z(-), 
the construction of the process R(-) in section 2 shows that an analogous change of 
variables formula holds for the process R(-) . Applying this formula to the right-hand 
side of (4.5) and differentiating under the integral, we obtain from the martingale 
property of the process in ( |4.5[ ) the heat equation 



(4.6) 



j n n 

d t p(t, r,r) = -w^vt d 2 Tk p(t, r, f) + ^ b k d Tk p(t, r, f) 



k=l 



k=l 



and the elastic boundary condition 

(4.7) q k d rk pit, r, r) - d rk+1 p(t, r, f) 

for t > 0, (r,f) G (W n ) 2 . 



whenever r k = r k+ \ 



Substituting the expression of (4.1) into (4.7), we deduce 



- ii e ^ n f iMi \^ a{l) -T t )D 2 f k+i ^ k+i \t^ a[k+1) -r k+ ^ = o 

<re§„ l^k+l 

whenever r k = r k+ i . Here, D 2 denotes differentiation with respect to the second 
argument. Plugging in r k for r k+ i in this last expression, and grouping together 
functions that have the same arguments, results in 

(4.8) 

Vl<Kn,l<4^2<i, tr(k) = l u a{k + 1) = £ 2 , a(k) = £ 2 , a(k + 1) = i x : 

= qt Ktr D 2 f k *(t,r ei -r k )f k+1 > e *(t,r h 



r k ) 



+ q k Kzf k+1 ^(t,r ei -r k )D 2 f k ^(t,r e2 -r k ) 
-qt +1 ^f kA (t,r h -r k )D 2 f k+1 ^(t,r £2 -r k ) 
- qt +1 ^D 2 f k+1 ^(t,r k - r k ) f k > e *(t,r h - r k ). 



Let us recall now (1.2), and take the Fourier transform with respect to the variables 
Tt x ~ r k (parameter a) and fg 2 — r k (parameter b), to obtain equations of the form 



(4.9) 

or equivalently 
(4.10) 



q (a G(t, a) H(t, b) + K(t, a) b L(t, b)) 

= (1-q) (G(t, a) b H(t, b)+a K(t, a) L(t, b)) 



G(t,a)H(t,b) 



qa + (1 - q) b 
q b + (1 — q) a 



K(t,a)L(t,b) 
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The only cases in which the fraction 

qa + (1 - q) b 
qb+ (1 — q) a 

can be factored as the product of a function only of a and a function only of b, are 
given by q = 1/2, q — 1 and q = 0. Moreover, up to multiplicative constants, 

(4.11) fk,e = fk+i,i^ e=l,2,...,n whenever q k = 1/2 , 

(4.12) D 2 f k/ = f k+1 ' e , £ = l,2,...,n whenever q~ = 1 , 

(4.13) f k ' £ = D 2 f k+l ' £ , £ = l,2,...,n whenever q~ = . 

Next, we introduce the linear parabolic operators 

(4.14) TZ k = d t -^a 2 k d 2 rk -b k d rk , k = l,2,...,n 



and substitute the expression of ( |4.1[ ) into (4.6), to obtain 

n 

(4.15) = ^ E (II -r<)) ftfc/^ftM*) -r») 

creS„ fc=l ^A; 

which shows 

(4.16) ft fc /*> CT W = 0, fc = l,2,...,n, <7G§„. 



This and (4.11 )-(4.13) imply b\ — b 2 — ■ ■ ■ — b n and (j\ — a 2 — ■ ■ ■ — a n . This 
completes the proof of the "only if" part. 

Conversely, if the conditions on the drift, dispersion and collision parameters in 
the proposition are satisfied, we define f k,i , 1 < k, £ < n by 

(4.17) f k ' k = cp bl;ai , fc = l,2,...,n, 

(4.18) f,t = f+i,e^ £=l,2,...,n whenever q k = 1/2 , k = 1, 2, n - 1 , 

(4.19) D 2 f k ' e = f k+1/ , £=l,2,...,n whenever q~ = 1, jfe = 1, 2, n - 1 , 

(4.20) f k >* = D 2 f k+1 ' e , £=l,2,...,n whenever = 0, jfe = 1, 2, n - 1 . 

We can express this state of affairs as follows: In order to determine the entry 
f k ' e (t, ■ ) in (4.1 ) for £ > k, we count the number u of ones in {q k , . . . , qj_{\ and the 
number 3 of zeros in {q k , . . . , qj_{\; and then compute f k ' e (t, ■ ) by differentiating the 
Gaussian probability density function <£& 1)01 (t, •) with respect to its second coordinate 
u times, and integrating the result with respect to the second coordinate 3 times. 
The entries f ' (t, ■ ) for £ < k are computed similarly. 
Now, in the case 

(4.21) <fc=V2, fc = l,2,...,n, 

one just needs to argue as in the beginning of this section to finish the proof; whereas 
in all other cases one can argue as in the proofs of Proposition 8 and Lemma 7 in [25] 
to deduce that the expression of (4.1) with these choices of f k ' e (t, •) , 1 < k, £ < n 
and K a = 1, it G S n gives the transition densities of the process R(-) . This completes 
the proof. □ 
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